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1. Introduction 

Accurate modeling of the mechanical and thermal properties of very-low- 
mass stars (VLMS), Brown Dwarfs (BD) and Extra-solar Giant Planets 
(EGP) is of prior importance for a wide range of physical and astrophysical 
problems, from the fundamental physics point of view to the astrophysical 
and cosmological implications. They provide natural laboratories to test 
the different theories, equations of state, nuclear reaction rates, model at- 
mospheres aimed at describing the physics of dense and cool objects. They 
represent the largest stellar population in the Galaxy, and thus provide a 
substantial contribution to the Galactic (disk) mass budget. Finally they 
represent one of the most intriguing questions in our understanding of the 
formation of star-like objects: are planet and star formation processes really 
different ? Is there, and if so what is, a minimum mass for the formation 
of star-like objects ? This field has blossomed recently with the discovery 
of several brown dwarfs (Nakajima et al. 1995; Rebolo et al., 1995) and 
numerous exoplanets since 51 Pegasi (Mayor and Queloz 1995; Mayor, this 
conference), which provide important information to challenge the theory. 

2. Theoretical improvements 

VLMS, or M-dwarfs, defined hereafter as objects with masses below O.6M , 
are compact objects, with characteristic radii in the range ~ 0.1 — 0.6 Rq. 
Their central densities and temperatures are respectively of the order of 
p c « 10 — 10 3 gem -3 and T c « 10 6 — 10 7 K, so that the stellar interior 
represents a strongly correlated plasma. BD's are defined as objects not 
massive enough for their central temperature to sustain hydrogen burning. 
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This characteristic stems from the onset of degeneracy in the contract- 
ing protostar, which prevents an increase of the temperature as contrac- 
tion proceeds. A BD therefore never reaches thermal equilibrium and cools 
down for its whole life. The hydrogen-burning minimum mass (HBMM) 
is ~ 0.07 — O.O9M , depending on the initial composition from Z=Z to 
Z=0, respectively (Chabrier & Baraffe, 1997). The minimum mass depends 
on stellar formation theory and is still uncertain (Hubbard, 1994). Brown 
dwarfs bridge the gap in the observed mass-distribution of astrophysical 
objects between the lowest-mass star and the largest solar planet, Jupiter 
(Mj = 0.001 M ). The distinction between BD and Giant Planets is based 
on different formation scenarios, but their inner structure (apart from the 
central rocky core in planets) and spectral signatures are both governed by 
the same physics. 

The effective temperatures of these low-mass objects are below T e ff ~ 
5000 K, down to ~ 100 K for BDs and EGPs, and surface gravities g = 
GM+/R% are in the range log g « 3.5 — 5.5. The low effective temperature 
allows the presence of stable molecules in the atmosphere (see Baraffe & 
Allard, these proceedings). The presence of these bands complicates tremen- 
dously the treatment of radiative transfer, not only because of the numerous 
transitions to include in the calculations, but also because the molecular ab- 
sorption coefficients strongly depend on the frequency. Moreover molecular 
recombination in the interior (2H — > H2) leads to a decrease of the adiabatic 
gradient so that convection penetrates deeply into the optically-thin atmo- 
spheric layers. Because of these physical processes, the grey-approximation 
is no longer valid below T e ff ~ 5000 K, as shown by several authors (Al- 
lard, 1995; Saumon et al., 1994; Chabrier et al., 1996; Chabrier &; Baraffe, 
1997). 

These conditions show that the modeling of low-mass objects requires 
a correct description of non-ideal effects in the interior, as well for the 
equation of state (EOS) as for the screening factors of the nuclear rates, 
a derivation of accurate models for dense and cool atmospheres, where 
molecular opacity becomes eventually the main source of absorption, and 
consistent (non-grey) boundary conditions between the atmosphere and the 
interior along evolution. 

A new EOS, devoted to the description of low- mass stars and giant 
planets, has been derived recently, which presents a consistent treatment of 
pressure ionization (Saumon, Chabrier & VanHorn 1995; SCVH), altough 
based on the so-called ideal-volume mixing between hydrogen and helium 
(see SCVH). Improvement in the field of cool atmosphere models has blos- 
somed within the past few years with the work of Allard and Hauschildt 
(1995; 1997), Brett (1995), and Saumon et al. (1994). Non-grey atmosphere 
models now extend down to T e jf = 900 K, since the discovery of the first 
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cool BD GL 299B offered a stringent test for such models (Allard et al. 
1996; Tsuji et al. 1996a; Marley et ah, 1996). 

A detailed description of the recent improvement in the theory of low- 
mass objects (LMS and BDs), and a comparison with standard grey-like 
treatments, is given in a recent paper by Chabrier & Baraffe (1997). The 
modelization of VLMS, BDs and EGPs and the confrontation with obser- 
vations is examined below. 

2.1. VLM STARS 

In spite of considerable progress in stellar theory - internal structure, model 
atmospheres and evolution - all the VLMS models so far failed to reproduce 
accurately the observed color-magnitude diagrams (CMD) of disk or halo 
stars below ~ 4000 K, i.e. ~ 0.4 — 0.6 M Q , depending on the metallicity. All 
the models predicted too hot an effective temperature for a given luminos- 
ity, i.e. were too blue compared to the observations by at least one mag- 
nitude. Such a disagreement stemmed essentially from shortcomings both 
in the physics of the interior, i.e. equation of state (EOS) and thus mass- 
radius relationship and adiabatic gradient, and in the atmosphere, since all 
models were based on grey atmospheres and approximate outer boundary 
conditions. Important progress has been made recently in this field with 
the derivation of evolutionary models based on a consistent treatment be- 
tween the interior and the atmosphere profile (Baraffe, Chabrier, Allard & 
Hauschildt 1995, 1997; Chabrier, Baraffe k Plez 1996). The effect of the 
outer boundary condition on the mass-T e jj relationship has been examined 
in detail by Chabrier & Baraffe (1997). These authors have made compar- 
ison with models based on different grey-like treatments and have shown 
convincingly that such treatments (which imply a T(r) relationship as an 
external boundary condition) are incorrect, or at best highly unreliable as 
soon as molecular formation sets in, i.e. for any object below T e ff ^ 5000 
K, i.e. m ~ 0.5 M . 

The LMS models based on the afore-mentioned updated physics and 
consistent (non-grey) boundary condition now reach quantitative agree- 
ment with observations for both the disk and the halo stellar population 
down to the bottom of the main sequence (Figure 1). For the disk popula- 
tion the models reproduce the observed color-magnitude diagrams both in 
the infrared (Mk vs (I-K)) and in the optical (My vs (V-I)) (Baraffe et ah, 
1997b), although below 0.1 M Q , My <; 15, the models are still too blue by 
~ 0.2 — 0.4 mag in (V-I), as illustrated in Figure 1. This is likely to stem 
from the still inaccurate TiO line list which shapes the flux in the opti- 
cal, or possibly from the onset of grain formation for solar metallicity (see 
Baraffe k, Allard, these proceedings). More importantly, the models are in 
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Figure 1. Mv-(V-I) diagram for different metallicities : [M/H]=-f.5 (solid line), 
[M/H]=-0.5 (dash-dot) and [M/H]=0 (dash). The Main Sequence of the Globular Cluster 
NGC6397 from Cool et al. (1996) is shown on the left part. Subdwarf halo field stars 
from Monet et al. (1992) arc indicated by full circles, as well as disk M-dwarfs of Monet 
et al. (1992) (crosses) and Dahn et al. (1995) (dots on the right hand side) 



excellent agreement with the observationally- determined mass-magnitude 
relationship (Henry & McCarthy, 1993) both in the infrared and in the op- 
tical (Chabrier, Baraffe & Plez, 1996; Baraffe et al., 1997b) as shown in the 
present Figure 2 and in Figure 1 of Baraffe & Allard (these proceedings). 
For metal- depleted populations, the models are in remarkable agreement 
with the main sequences of globular clusters observed with the HST nearly 
down to the bottom of the main sequence and with the halo field subdwarf 
sequence, as shown in Figure 1 (Baraffe et al., 1997a). We stress that, since 
LMS are essentially fully convective below ~ 0.4 M & (Chabrier & Baraffe, 
1997), the models are not hampered by any adjustable parameter and the 
agreement between theory and observation reflects directly the reliability 
of the physics entering the models. These models yield the derivation of 
reliable mass-functions for the disk population (Mera, Chabrier & Baraffe, 
1996) and for globular clusters and halo field stars down to the brown dwarf 
limit (Chabrier and Mera, 1997). 
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Figure 2. Mass-luminosity relationship in the K-baxiA. The dash-line is the fit of the 
observationally-determined relation (Henry & McCarthy, 1993; HMc93), the solid line is 
based on the most recent models for solar metallicity (Baraffe et al., 1997b) for t = 5 
Gyr, whereas the dash-dot line is the same models for t — 0.1 Gyr (the ZAMS age for 
0.1 Mq is about 0.5 Gyr). Note that the three lowest-mass objects observed by HMc93, 
shown on the figure with their error bars, all exhibit strong surface (H a -emission) and 
coronal (X-emission) activity. Although this is not a proof of a young age, it is a strong 
hint and is supported by the present analysis. 

2.2. BROWN DWARFS 

The recent discovery of substellar objects enables us to test the reliability of 
the theory in the BD regime. The first young brown dwarfs were identified 
in the Pleiades star cluster (Rebolo et al. 1995), namely Teide 1 and Calar 
3, which both have retained their initial Lithium abundance (Rebolo et 
al., 1996). Their extremely low luminosity (log L/Lq ~ —3.1) combined 
with the presence of lithium yields masses m ~ O.O5M , well below the 
hydrogen-burning minimum mass (Rebolo et al., 1996; Baraffe et al, 1997b). 
Given the large effective temperature of these two BDs {T e ff ~ 2500K), 
due to their young age (~ 10 s yrs), the most stringent test for extremely 
cool models is provided by Gl 229B (Nakajima et al. 1996). The presence 
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of methane in its infrared spectrum clearly reveals its substellar nature, 
since CH4 absorption appears below 1600 K (Tsuji et al., 1996a) and a 
star at the hydrogen burning limit has an effective temperature of about 
2000 K (Baraffe et al. 1995). Comparison of synthetic and observed spectra 
yields an effective temperature of ~ 1000-K" (Allard et al. 1996; Tsuji et 
al. 1996a; Marley et al. 1996). Uncertainties in the age of the system and 
in the temperature of G1229B yield some indetermination for the mass. 
Evolutionary calculations based on the afore-mentioned synthetic spectra 
and non-grey atmosphere models yield the most likely solution M ss 0.04 — 
0.055 M Q for an age ~ 5 Gyr, similar to our solar system (Allard et al. 
1996; Marley et al., 1996). 

Figure 3 presents different isochrones as a function of metallicity (Baraffe 
et al. 1997b) at the bottom and below the main sequence in IR colors. The 
K-limit magnitudes corresponding to the HBMM are indicated by full cir- 
cles. The blue loop displayed in the substellar domain stems from collision- 
induced absorption of H2 , and CH4 absorption for solar-metallicity, in the 
K-band (Saumon et al., 1994; Baraffe et al., 1997a,b) which shifts the flux 
back to shorter wavelengths. We predict this blueshift in IR-colors, whereas 
optical colors keep reddening almost linearly, to be the most important pho- 
tometric signature of the transition from the stellar to the sub-stellar do- 
main. We have made predictions of such signatures in the NICMOS filters, 
which should be verifiable in a very-near future (Baraffe et al., 1997a). 

2.3. EXTRASOLAR GIANT PLANETS 

A general theory of EGPs, from 0.3 to 15 Jupiter masses, and the possi- 
bility of detection by existing and future observational projects, has been 
performed recently by Saumon et al. (1996). These calculations take root in 
the theory of solar giant planets derived previously (Chabrier et al., 1992; 
Guillot et al. 1995). In this first generation of calculations, the spectral 
emission of the EGP and of the parent star is approximated by a black 
body distribution. Although weakened by this approximation, the predic- 
tions made in Saumon et al. (1996) represent the first benchmark in search 
strategies for the detection of EGPs. They have been applied more specif- 
ically to the case of 51 PegA to demonstrate the stability of giant planets 
orbiting nearby stars (Guillot et al., 1996). 

More recent calculations by Marley et al. and Allard et al. do now 
predict synthetic spectra in the substellar domain from T e ff ~ 2000 K 
down to ~ 300 K. The predicted absolute fluxes of BD or EGP (Allard et 
al. 1997) are displayed in Fig. 3 of Baraffe & Allard (these proceedings) and 
compared to the sensitivity of ground and space-based observing platforms. 

An extra-degree of complication in the atmosphere of objects below 
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Figure 3. Color-magnitude diagram in the near IR for the metallicity and ages indicated 
on the figure. The full circles on the curve indicate the stellar/sub-stellar transition. 



T e ff ~ 2600i^ is due to the onset of grain formation, as suggested by Tsuji 
and collaborators (Tsuji et al. 1996b). This is illustrated in Figure 4 which 
displays the effective temperature and the radius of various exoplanets dis- 
covered recently with the domain of condensation (bars on the right hand 
side) of various compounds (Guillot et ah, 1996). 

3. Conclusion 

The recent improvements in the description of the mechanical and thermal 
properties of cool, compact objects, and of their photometric signature 
now provide solid grounds to analyse the observations and make reliable 
predictions. The three essential inputs in the theory of these objects are an 
accurate EOS, with a reliable treatment of non-ideal effects and pressure 
ionization, synthetic spectra with accurate molecular absorption coefficients 
in the optical and the IR, and consistent evolutionary calculations with 
correct (non-grey) boundary conditions between the atmosphere and the 
interior profiles. 

This improved theory of LMS and substellar objects allows now the 
derivation of reliable stellar mass-functions down to the brown dwarf limit 
and, depending on the rate of discovery of substellar objects in the near 
future, will yield eventually the brown dwarf mass function, in connection 
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Figure 4- T e ff and radius as a function of mass for several EGPs (kindly provided by 
T. Guillot & A. Burrows). Jupiter and Saturn are also indicated by their usual symbols. 
The bars on the right hand side show the domain of condensation of different species. 



with microlensing observations. This will bring new insight on the still 
unsolved problem of stellar versus planetary formation. 

The derivation of cool atmosphere models and synthetic spectra includ- 
ing grain formation represents the next substantial improvement. Work is 
already under progress in this field, both in the Tucson group and in the 
Lyon- Wichita group, and promises more exciting results in the near future. 
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DISCUSSION 

PIERRE MAXTED: You suggested that two of the brown dwarfs seen in 
the Pleiades are binary stars since they sit so far above the 
isochrones in the HR diagram. Is there any independent 
observational evidence of binarity in these stars? 

GILLES CHABRIER: The groups of Rafael Rebolo and Gibor Basri, who 
suggested the possible binarity, are examinating this possibility 
but, as far as I know, they don't have a definitive answer yet. 

J0RGEN CHRISTENSEN-DALSGAARD : Would it be fair to ask what the 
MACHO objects are, if not brown dwarfs? 

GILLES CHABRIER: As we've shown in a recent paper (ApJ 468 L21) , 
they can be explained by halo white dwarfs, providing some 
stringent conditions are met. I will illustrate this point in my 
white dwarf talk. 

FLAVI0 FUSI PECCI: You have shown nice fits of observed faint main 
sequences with Galactic globulars and claim that the quality of the 
fits is much better than obtained so far, thanks to the improved 
quality of your models compared to previous ones. I have seen 
similar fits obtained, for instance, for NGC 6397 (HST Data) by 
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Alexander et al . (1997 A&A, 317 90) using their own models - could 
you schematically explain which are the main differences between 
the two sets and how you determine and evaluate the quality of the 
best fitting solution? 

GILLES CHABRIER: A detailed answer is given in our papers (Chabrier 
& Baraffe 1997; Baraffe et al . 1997a). To summarize, there are five 
important differences: 

1. Our models rely on consistent non-grey boundary conditions, 
whereas the Teramo models use a grey condition, based on T(t) 
approximation. Even when modified (or adjusted), the grey 
approximation is basically not correct/reliable for LMS, 
because of the very physics characteristic of these objects 
(see the present review) . 

2. With the same physics, and no adjustable parameter, our models 
reproduce all the observed sequences from Fe/H = -2.2 (M 15) to 
-1.0 (lu Cen) and even -0.5 (47 Tuc) , not just NGC 6397. 

3. The Teramo models assume [M/H] = [Fe/H] which is not correct 
for metal-depleted stars, because of the O/Fe enrichment, 
which yields [M/H] = [O/H] . 

4. The Teramo models rely on previous (so-called "Base") colors by 
Allard & Hauschildt which overestimate the opacity (because of 
the Straight mean approximation, see Allard et al . , 1997, 
ARA&A) and are now abandoned. Moreover the bolometric 
corrections have been shifted arbitrarily to recover the Kurucz 
values at high temperatures (see their §4) . 

5. Our models use consistent reddening corrections, calculated 
from Allard' s synthetic spectra. These corrections 
(E(V-I)=0 . 22) are in excellent agreement with the ones used by 
the observers (E(V-I)=0 . 23; Cool et al . , 1996). The ones used 
by the Teramo group (E(V-I)=0. 19, see their Fig. 7) differ 
significantly. With the reddening and metallicity quoted by the 
observers, their models are too red w.r.t. the observations 
(see e.g. Fig 5 of Baraffe et al . , 1997a). 



